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The Starobinsky model is a natural inflationary scenario in which inflation arises due 
to quantum effects of the massless matter fields. A modified version of the Starobinsky 
(MSt) model takes the masses of matter fields and the cosmological constant, A, into 
account. The equations of motion become much more complicated however approximate 
analytic and numeric solutions are possible. In the MSt model, inflation starts due to the 
supersymmetric (SUSY) particle content of the underlying theory and the transition to 
the radiation dominated epoch occurs due to the relatively heavy s-particles decoupling. 
For A = the inflationary solution is stable until the last stage, just before decoupling. 
In the present paper we generalize this result for A ^ 0, since A should be non-vanishing 
at the SUSY scale. We also take into account the radiative corrections to A. The main 
result is that the inflationary solution of the MSt model remains robust and stable. 

Keywords: effective action; inflation; conformal anomaly; supersymmetry. 



1. Introduction 

The standard cosmological model covers a wide class of phenomena and fits the 
c urren t observational tests with great success. However, this model has problems 
the initial singularity, horizon, flatness and monopoles in the early period 
of the universe. These problems can be solved if we assume that the primordial 
universe starts with a very fast expansion, denominated inflation by Guth in 1981 

An essential natural inflationary scenario is one in which inflation is driven by 
quantum corrections to the Einstein-Hilbert action, suggested by Starobinsky in 
1980.^ The Starobinsky model is based on the semiclassical approach to quantum 
field theory (QFT) in curved space-time. Within this theory the metric is treated as 
a classical background for the quantum dynamics of the matter fields. This approach 
presents a c onsi stent theory at energies of a few orders of magnitude below the 
Planck scaled 

In the original Starobinsky model, inflation is a consequence of the quantum 
effects of massless matter fields.™ The model assumes a non-minimal conformal 
coupling between the scalar field and gravity, £ = 1/6. In this case, the mass- 
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less matter fields are conformally invariant having a traceless stress tensor at the 
classical level. However, the one-loop contributions create a trace anomaly which 
changes the dynamics of the conformal f actor of the metric (see Refs. [4j [TJ [8]) and 
also the metric and density perturbations E PO I H I An alternative option is to apply 
the effective action method, using the conformal anomaly to calculate the induced 
effective action.™! Infkt ion naturally arises from the total action which is ob- 
tained from the sum of the anomaly-induced effective action to the classical terms, 



including the Einstein-Hilbert one!- 



3 pIHT8] 

An alternative version of the Starobinsky model was proposed in Refs. Qj|l [Ml 
[2T1 The main advantage of this modified version is that inflation starts in the sta- 
ble regi me whi ch is afterwords interpolating to an unstable regime at the end of 
inflationSEELZlTlie modified Starobinsky (MSt) model is a natural extension of the 
Starobinsky model. In the MSt version, inflation is due to the contribution of the 
quantum effects of both massless conformal and massive matter fieldsP^EH xhe 
massive theory is not conformally invariant at the classical level due to the masses 
of the scalar and fermion fields. However, using a conformal description, the massive 
matter fields become conformally invariant and we can use the conformal anomaly 
method to derive the effective actionf^QEIEl] 

The stability condition depends on the particle content of the underlying quan- 
tum field theory. Assuming supersymmetry in the high energy region, the super- 
symmetric (SUSY) particle content provides an initial inflationary period stable 
under small perturbations of the conformal factor. 1 ^ The Hubble parameter, H, is 
not constant at this stage as in the original Starobinsky model. Instead, the infla- 
tionary expansion is slowing down due to the contributions of massive particles. At 
some point the stable inflation becomes unstable when the s-particles decouple and 
supersymmetry breaks down. 

For a vanishing cosmological constant (CC) we showed that the numerical solu- 
tion of the MSt model has an accur ate app roximate analytic solution which is robust 
during the entire inflationary stageP2l21] Recent astronomical observations indicate 
a small value for the CC at present, A present « 0.7/3° ~ 10 -121 Mpy . There is a large 
discrepancy comparing A prcsont with the vacuum energy density p V ac = A/(8ttG) 
at higher energy scales, where Mpi = 1/V8ttG = 2.44 x 10 18 GeV is the reduced 
Planck mass. This discrepancy is known nowadays as the old cosmological problem 
(see Ref . [23l for a classical review). In fact, p V ac ~ ^susy — 10~ 12 Mpj for example, 
on the SUSY scale. The value of p v&c is due to the extremely exact fine-tuning of 
the vacuum counterpart of the CC today and to the abrupt change of the CC due 
to the induced counterpart, which presumably took place at an early stage of the 
evolution of the universe. One can find a discussion of the CC problem in the QFT 
framework in Ref. 1291 

In the MSt model, the contribution of the massive scalar fields emerges in the 
effective action of gravity through the renormalization of the cosmological constant 
A term in Einstein's equations. Hence, when we take A = in Eq.Q, the MSt 
model is left solely with the contribution due to the masses of the fermion fields, 
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which is within the definition of / [Eq. (|10p ]. the /3-function which renormalizes the 
gravitational constant. 

On the other hand, the /3-function which renormalizes the A term can be linked 
to a dimensionless expression g (see Eq. (fTTj) below). The parameter g is given by 
an algebraic sum of the fourth powers of the masses of the fermion and scalar fields, 
taking their statistic and multiplicities into account. It is not obvious that this sum 
should cancel out at all stages of inflation, even if the supersymmetry is initially 
present. Therefore g may not vanish. Hence, when we take A ^ the contribution of 
the massive scalar fields should, in principle, contribute to the solution of the MSt 
model and the effect of such contributions on the inflation should be investigated. 

The stability criterion in the MSt model depends, in principle, as we are going 
to show, on the size and sign of g and on the supersymmetry breaking (the value 
of A). The stability at the initial inflationary period is one of the great successes 
of the MSt model. The MSt model does not present any initial condition problem 
for A = since the stability criterion is satisfiedPH In this paper we reconsider 
the solutions and the stability condition for the MSt model assuming the minimum 
supersymmetric standard model (MSSM) particle content at the beginning of infla- 
tion and the natural values of g and A at the corresponding energy scale. We show 
here that, in contrast to the naive expectations, a non- vanishing CC and g (the last 
depending of multiplicities and masses of the fermion and scalar constituents of the 
supersymmetric model) do not destroy the stability in the MSt model. 

The paper is organized as follows. In section [2l we present the framework of 
the modified Starobinsky (MSt) model, revisiting the stability criterion for a non- 
vanishing A and g as well as the approximate analytic solution for the MSt model in 
section[3l We introduce the natural values of the parameters to analyze the solutions 
and stability criterion numerically in section 21 The conclusions of the paper are 
presented in section [5l 

2. The framework of the modified Starobinsky model 

In this section, we discuss the anomaly induced inflation formalism of the Starobin- 
sky models following the notations in Refs. [T?ll2T)l In the Starobinsky model, infla- 
tion comes from the contr ibution of the quantum effects of massless matter (scalar, 
fermion and vector) fields SH3IS] Assuming a non-minimal conformal coupling be- 
tween the scalar field and gravity, £ = 1/6, the massless matter fields are invariant 
under the local conformal transformation of the fields and the metric 

g^^g^e 2 ^, (1) 

where cr(rf) = In a(rj), dt = adrj, rj is the proper (conformal) time and a is the scale 
factor. 

The massless matter action satisfies the conformal Noether identity which im- 
plies that the energy-momentum tensor is traceless at the classical level. How- 
ever, the one-loop quantum contribution of the massless matter fields cause a trace 
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anomaly < T >=/= which is useful to calculate the induced effective actionfSH^ISl 
In this conformal anomaly method, the anomaly-induced effective action Ti nc i is 
derived from the trace anomaly using the conformal metric ([T]). We assume the 
spatially flat (k = 0) line element of the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) metric 

ds 2 = a 2 (r])(dT] 2 -dx 2 ). (2) 

General equations for k = ±1 (non flat space) can be found in Refs. HI [IS] for the 
massless theory and for the modified Starobinsky (MSt) model in Ref. [2TJ 

In the MSt model, inflation starts with a massive supersymmetric (SUSY) par- 
ticle content iVr^i/2,1, where 0, 1/2 and 1 correspond to scalar, fermion and vector 
fields, respectively. The theory is not conformal invariant anymore due to the ma sses 
of the scalar and fermion fields. However, applying a conformal description}^ the 
massive theory becomes conformal invariant at the classical level and we can use 
the conformal anomaly method to derive the effective action.^EB 

The classical vacuum action which provides the possibility to renormalize the 
massive theory, in its minimal version contains, besides the Newton G and A terms, 
the parameters a\, 02, 0.3, defined according to 

Svac = ShD + SeH , 

where Shd is the part which contains higher derivatives of the metric@ 

S HD = J d A x^~g { ai C 2 + a 2 E + a 3 V 2 i?} , (3) 

where C 2 = C^^C^^ = Rp Vaf} W lva f ) - 2R aP R afi + ± R 2 and E = 
R^vapR^"^ — 4 R a fiR al3 + R 2 are the square of the Weyl tensor and the integrand 
of the Gauss-Bonnet topological term, respectively, and Seh is the Einstein-Hilbert 
action 

1 r ^ 



Seh = - ^ / d^xj—g (R + 2A) . (4) 

Using the method proposed in Refs. [2111 |2"T1 the anomaly-induced effective action 
Find can be derived from the trace anomaly of the massive theory with the usual 
conformal transformations to the fields and the conformal metric (QJ. We obtain 
the total effective action adding the classical terms 

T = Shd + J d 4 xy/^{wC 2 a + b(E-^V 2 R)a + 2baAa} 



3c + 26 
36 



d 4 x^g~R 2 - / d A x^e Aa 



A 
8^G 



9 ■ o 



a One has to notice that the introduction of the non-conformal terms like f y/—gR 2 is possible, 
but not necessary, for the renormalization of the free conformal invariant theories (see e.g. Ref ll7l 
and references therein). 
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- / d A X, 



-ge 2 ° [R + 6(Va) 2 } 



1 



16ttG 



(5) 



where we have discarded a possible unknown term which comes from the integration 
of ^ind- The coefficients w, b, c, / and g are the /3-functions for the parameters of 
the classical vacuum action ai, 03,03, G and A, respectively. The explicit form of 
these functions is: 



l 
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2 



(4tt) 2 



(6) 

(7) 
(8) 



Here Nf and iVs are the multiplicities of the fermion and scalar massive fields with 
the masses rrif and ms- Taking the minimal variation of Eq.© with respect to the 
conformal factor a, we obtain the following equation of motion 

"a (t) + 7ct (t)<r (t) + 4 cr 2 (i) + 4 f 3 - ' 



(a (t) + 2a 2 (t) 



a(t)& (t) 
Mt) 



c 



\f* (*) 



+ 



2 AM|j 
3c 



1 - g(7(t) 



0. 



(9) 



where M P1 = l/y/8nG = 2.44 x 10 18 GeV is the reduced Planck mass. The param- 
eters / and g are defined as dimensionlcss functions of the previous / and g 

/ = (16ttG) / 



1 



9 



A/(8ttG) 2(47t) : 



E 



3(4tt) ,. 
N s m 4 s 



N f m 2 



(10) 



pi 



(47T) 



E 



f 



A ad M 



(11) 



PI 



where we have introduced here a dimensionless vacuum parameter, according to 
Pvac = A/(8ttG), as 

A ad = /w/MpV (12) 

The Starobinsky inflationary solution can be found assuming the massless matter 
fields (/ = g = 0) and A = in Eq.© 



a(t) 



„o-s(*) 



(13) 
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where Hs = H a Mpi is the Hubble parameter in the original Starobinsky models 
and H a = 1/^—2 b is a dimensionless constant. The parameter H a depends on the 
particle content A , Ni/ 2 and Ni according to Eq.(j6|). A solution for A^O can be 
found in Ref. [H 

When massive field contributions are taken into account, the equation of mo- 
tion (|9|) becomes more complicated and can not be solved analytically. N umeric al 
calculations with A = have shown that H is slowly decreasing in timePSEH ft 
is important that the corresponding solution is stable until the transition point, 
where heavy s-particles decouple, supersymmetry breaks down and another phase 
of inflation startsP During the stable period one can find an accurate approximate 
solution which closely reproduces the numerical solution. In the next two sections 
we shall generalize these results by taking the cosmological constant and its running 
into account. 



3. The stability conditions and the approximate solution 

In the massless theory, the stability condition is well-known: c > 0.SE2 This condi- 
tion corresponds to the assymptotic stability of the de Sitter solution under small 
perturbations of the conformal factor of the metric cr(r) — > cts(t) + y(r), where as 
is given by Eq. (|T3"]) . According to Eq.flSJ), this condition corresponds to the following 
inequality for the particle content of the theory 



1 



1 



Ai < -A 1/2 + — A , 



(14) 



where Ai, N\n, Nq are the numbers of particles with the corresponding spin. 

Furthermore, the stability condition for the massive theory, which was obtained 
in Ref. 1211 has the following form 



3/ff 2 (r)-5A(r) > 0, 



with 



H\r) = - 



(r) 



4b 



1 



8b A(? 



1 + ^r — 



3 M|,(r) 



1/2 



(15) 



(16) 



where we have normalized the time to r 
time) and 

MUt) = Ml, \l - fa(r) 



t/tpi (tpi ~ 5.3 x 10 _44 sec is the Planck 



A(r)=A{l-5a(r)~g/4} 



(17) 



During inflation, the solution of the MSt model is stable (obeying the criterion of 
stability) until a characteristic scale M* After this scale, or the dimensionless 
scale n = M^/Mpi, the Hubble parameter decreases, becoming constant and very 
small = (t(t) w /i, when the sparticles decouple and the matter content A 0!l / 21 
becomes modified. As a result, the inequality in Eq. (ll4p changes sign and the uni- 
verse enters into an unstable inflation regime with an eventual transition to the 
FRW evolution. We have shown that the last transition of the inflation, satisfying 
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the stability condition (fT4"]) , is independent of the value of the cosmological constant 
A or the curvature k in Ref. [211 

Our purpose here is to consider the MSt model assuming the possible contribu- 
tion of the cosmological constant A and g to the numerical solution. The approx- 
imate solution for the MSt model is obt ained a ssuming / small, a slowly varying 
Hubble parameter and, in addition, A = O .HQED] Discarding the term / g 2 in Eq .([9l) 
we obtain exactly the equation of motion for the massless matter fieldsP^D^ The 
approximate analytic solution is then obtained taking A = in Eq. (|17p . Thus, ac- 
cording to the first expression of EnJlT]), integrating the running Hubble parameter 
H(t) = cr'(r) = H a Mpi(t) we fincE 

a(r) = Ha r-^r\ (18) 

This approximate analytic solution can be understood as the massive contribution 
(with / 7^ and A = 0) arises as a running G, G(r) = [87tMpj(t)] -1 , to the 
Starobinsky model. The approximate solution Eq. (flBf shows a very good agreement 
with the numerical solution of the MSt model for the case A = OpU However, l e t us 
remark that the corresponding analysis was performed for the particular case A = 
and g = 0. In particular, we assumed A = and that our numerical solutions is 
stable as far as / > [Eq. (JTDJ)] , such that the criterion (fTSj) was satisfied. However, 
the stability criterion for the massive theory (| 1 5 j) may depend, in principle, on the 
value of A and also on the size and sign of g. The last parameter can be different 
from zero, since it is not necessary that there is an algebraic cancelation between 
the sum of the fourth powers of the masses of the fermion and scalar fields. Using 
a natural value for A during the inflationary period and the possible values of g, in 
the next section we analyze the solutions and the stability criterion for the general 
MSt model. 



4. Numerical analysis 

In this section, we discuss the natural values of the parameters / and g for A at 
the SUSY scale. Using these values we compare the approximate parabolic solution 
with the full numerical solution of the MSt model and the stability criterion for 
the initial evolution of the inflationary period. Recent astronomical observations 
indicate a small value for the cosmological constant 

A obs « 0.7p° c ~ 10" 47 GeV 4 ~ 10~ 121 . (19) 



b Similar behaviour arises in simple inflationary models with potential V <x </> 2 (see Ref. 131 1 and 
in the limiting form of the original Starobinsky model, called "quasi-de Sitter" stage in Ref. 
1111 Therefore, from the theoretical point of view, the advantage of the MSt model is that it is 
completely free of any ad hoc assumption (such as inflaton potential) and it starts with a very 
stable inflationary solution without any fine tuning of the parameters. In the massless theory, for 
instance, it is necessary to introduce some additional term or terms into the classical action of 
vacuum (e.g. a sufficiently large positive coefficient in the (J ^/~gR 2 )-term) to provide c > 0. 
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It can be much more larger at higher energy scales, for example, on SUSY (or 
GUT) scales. This difference can be due to phase transitions, e.g. the electroweak, 
which took place between the present time and the SUSY scale. Let us suppose 
Msusy ~ 10 16 GeV ~ 10~ 3 Af P1 which is consistent with the GUT scale (~ 10 14 - 
10 16 )GeV. In this case, 

Pvae-MlusY^lO- 12 ^/ 4 !. (20) 

The difference between the two above values is a manifestation of the well-known 
cosmological constant (CC) problem. Certainly it is much easier to assume that 
the vacuum energy is zero than to solve the CC problem.^2H2H] However, let us 
remark that almost all deviations from a "perfect equilibrium" vacuum lead to a 
non-zero vacuum energy density which is proportional to the perturbations of the 
vacuumjfl Nevertheless, from Eq. (|20|) , the natural value of the dimensionless vacuum 
parameter A a d [Eq.fTJ}] at the SUSY scale is 

A ad « 10~ 12 . (21) 

For the numerical analysis, we use the particle content iV^i/2,0 — (12,48,104) of 
the minimum supersymmetric standard model (MSSM) and / in the range 

/~(10- 5 -10- 3 ), (22) 

which corresponds to the scale of Eq. (|2T))) and a typical matter content Nf + N s ~ 
10 - 1000. 

Whether we assume the masses of the fermions and scalars to be of the same 
order or not, we can stipulate the parameter g from Eq. (|ll[) as 

Let us extrapolate the possible difference between the sum of the scalar masses 
squared and A a( j, as well as the simplified approximation used in Eq. (|23p . in the 
range (10~ 3 - 10 5 ). Thus, using Eqs.fl2I} and (22} we find 

g ~ ± (10~ 8 - 10 2 ) . (24) 

We begin the analysis solving the equation of motion of Eq. @ for the dimen- 
sionless A at the SUSY scale [Eq. (|2"Tj) ] varying / and g within the range of Eqs. (f2"2"|) 
and (f2~4")l . respectively. We check the stability criterion [Eq.(fr5}] during the infla- 
tionary period assuming the possible combinations of the range in / and g, as well 
as the possibility of positive and negative signs in g until the end of inflation which 
should happens at t cn d ~ 2/H a f, according to the approximate parabolic solution 
(fT8| . We checked that the stability criterion [Eq. (fT5| ] is satisfied for all possible 
combinations of the values of the parameters. We found that 3fH 2 (r) — 3 A(r) is 
positive until the end of inflation and the numerical solution is independent of the 



c For a recent discussion see e.g. Refs. 24 and 1251 
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Fig. 1. a) The stability criterion for the massive theory Eq, l|15|l assuming the MSSM particle 
content = (12,48, 104), the natural value of the dimensionless A at the SUSY scale (1 2 1 I t 

and a combination of extreme values of the parameters / = 10 — 3 and g = — 10 2 . b) The time 
dependence of A(r) from the transformations Eg. (1171 1 for the same values of the parameters as in 
Fig. la. 

initial conditions. As an example, we show the stability criterion for the massive 
theory assuming a MSSM particle content, the natural value of the dimensionless A 
at the SUSY scale [Eq.(|2ip] and a combination of extreme values of the parameters 
/ = 10 -3 and g = — 10 2 in Fig. la. The result can be understood looking at the 
time dependence of A(r) from the transformations of Eq.(|17p, We notice that one 
can neglect the constant terms so that the function A(t) has the opposite sign of g 
(since ct'(t) = H(r) > 0), or A(r) cx —gcr(T) as shown in Fig. lb. This mean that 
the stability criterion is then ~ 3fH 2 (r) +g 2 <j(t) which always has a positive sign. 

The numerical solution of Eq. (|9|) can be compared to the approximate parabolic 
solution [Eq.(|18[)] for all possible combinations of the parameters values discussed 
above. As a result we found an accurate compatibility between the numerical and 
approximate solutions. We show in Fig. 2 the approximate versus the numerical solu- 
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100 200 300 400 

T S t/t H 

Fig. 2. The numerical solution of Eq. for the same values of the parameters as in Fig.l versus 
the approximate parabolic solution (the dashed line). 

tion for the MSSM particle content with the same extreme values of the parameters 
used in Fig.l. The dashed line shows the approximate parabolic solution [Eq. (fT8)) ]. 
The contribution of the difference due to massive fermions and bosons, which is 
within the definition of g can not be distinguished in the numerical solutions. The 
massive fermions are mainly responsible for slowing down cr(r) and we can safely 
consider g and A vanishing as an approximation. The approximate solution is also 
robust in this stage. 

5. Conclusions 

We considered the inflationary solution in the modified Starobinsky (MSt) model, 
taking into account the cosmological constant A and quantum contributions to the 
vacuum energy density from massive fermions and bosons. The fields with differ- 
ent statistics contribute with opposite signs and the overall quantum effect is due 
to the difference of their contributions. The natural value for A corresponds to the 
supersymmetry breaking value at the MSSM scale. We also assumed the correspond- 
ing particle spectrum for evaluating the quantum contributions. The corresponding 
dimensionless parameter is called g. It turns out that the agreement between ap- 
proximate analytical and numerical solutions is robust for the natural values of the 
above-mentioned parameters. Furthermore, the numerical solution is not sensitive 
to the sign of g. The latter means that the contribution of the massive fermions can 
be smaller or larger than the contribution of the massive boson particles without 
visible effect to the numerical solution of the model. The inflation remains stable 
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until the point when the s-particles decouple and supersymmetry breaks down, for 
any natural choice of the parameters. 

On the other hand, the renormalization of the inverse Newton constant, 
parametrized by the dimensionless quantity /, is very important. This parameter 
is always positive and its magnitude depends only on the fermion spectrum of the 
theory. Indeed, / is responsible for decreasing the value of the Hubble parameter in 
the course of inflation. 

Finally, it is very important that the value of the cosmological constant A and 
the contributions of the massive particles do not destroy the stability which holds at 
the initial stage of inflation. Hence, the MSt model provides a possibility to describe 
all stages of inflation without fine-tuning of the parameters of the theory and/or 
initial data. 
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